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ABSTRACT 

A  comparison  was  made  between  computer  model  predictions  of 
gear  dynamic  behaviour  and  experimental  results.  The  experimental 
data  were  derived  from  the  NASA  gear  noise  rig,  which  was  used  to 
record  dynamic  tooth  loads  and  vibration.  The  experimental  results 
were  compared  with  predictions  from  the  Australian  Defence  Science 
and  Technology  Organisation  Aeronautical  Research  Laboratory’s 
gear  dynamics  code,  for  a  matrix  of  28  load-speed  points.  At  high 
torque  the  peak  dynamic  load  predictions  agree  with  experimental 
results  with  an  average  error  of  5  percent  in  the  speed  range  800  to 
6000  rpm.  Tooth  separation  (or  bounce),  which  was  observed  in  the 
experimental  data  for  light-torque,  high-speed  conditions,  was 
simulated  by  the  computer  model.  The  model  was  also  successful  in 
simulating  the  degree  of  load  sharing  between  gear  teeth  in  the 
multiple-tooth-contact  region. 

INTRODUCTION 

The  dynamic  behaviour  of  gears  is  important  for  a  number  of 
reasons — the  dynamic  load  increment,  the  gear  life,  the  gear  noise, 
and  the  overall  vibratory  behaviour  of  the  gear  system.  In  aero¬ 
space  applications,  such  as  helicopter  transmissions,  all  of  these 
factors  are  relevant,  and  weight  is  also  an  important  factor.  A 
knowledge  of  dynamic  load  factors  can  assist  in  weight  reduction. 

The  NASA  gear  noise  rig  was  built  to  enable  fundamental 
studies  of  gear  dynamic  behaviour  to  be  carried  out  and  to  provide 
support  to  gear  noise  reduction  programs.  Oswald  et  al.  (1991) 
compared  dynamic  load  measurements  from  the  NASA  rig  with  pre¬ 
dictions  from  the  computer  program  DANST  (Dynamic  Analysis  of 
Spur  Gear  Transmissions).  This  report  continues  that  work,  and 
compares  the  same  experimental  data  with  predictions  from  another 
model,  developed  at  the  Australian  Defence  Science  and  Technology 
Organisation  (DSTO)  Aeronautical  Research  Laboratory  (ARL), 
designated  ARL_DYN. 

The  two  models  differ  significantly.  DANST  (Lin  et  al.  1989, 
1987)  has  four  torsional  degrees  of  freedom  representing  the  motor, 
the  gears,  and  the  load.  The  shaft  stiffness  elements  are  represented 
by  linear  stiffnesses,  and  the  tooth  compliance  by  a  variable  stiff¬ 


ness.  ARL_DYN  (Rebbechi,  1991)  models  the  gear  in  a  way  that 
allows  the  detailed  contact  conditions  of  each  of  the  tooth  pairs  in 
contact  to  be  separately  represented.  The  model  also  includes  shaft 
deflection.  Although  this  has  the  disadvantage  of  increasing  the 
complexity  of  the  model,  it  has  the  advantage  that  the  model  para¬ 
meters,  such  as  the  sliding  friction  coefficient,  can  be  varied  in  a 
way  that  accords  with  the  physical  parameters. 

The  objectives  of  this  work  were  (1)  to  evaluate  predictions  of 
dynamic  load  'increment  by  using  the  ARL  model  and  compare 
these  with  experimental  results,  and  (2)  to  refine  the  model  as 
necessary  to  match  the  predictions  with  the  measured  dynamic  load 
curves. 


APPARATUS  AND  PROCEDURES 

The  NASA  Lewis  gear  noise  rig  (Fig.  l)  was  used  for  these 
tests.  This  rig  features  a  single-mesh  gearbox  powered  by  a  150-kW 
(200-hp)  variable-speed  electric  motor.  An  eddy-current  dynamo¬ 
meter  loads  the  output  shaft.  The  gearbox  can  be  operated  at 
speeds  up  to  6000  rpm.  The  rig  was  built  to  carry  out  fundamental 
studies  of  gear  noise  and  the  dynamic  behaviour  of  gear  systems.  It 
was  designed  to  allow  testing  of  various  configurations  of  gears, 
bearings,  dampers,  and  supports.  The  gearbox  is  extensively  instru¬ 
mented  for  strain  and  vibration  measurements. 

A  poiy-V  belt  drive  was  used  as  a  speed  increaser  between  the 
motor  and  the  input  shaft.  A  soft  coupling  was  installed  on  the 
input  shaft  to  reduce  input  torque  fluctuations,  which  were  caused 
by  nonuniformity  at  the  belt  splice. 

The  test  gears  were  identical  spur  gears  (at  1:1  ratio)  machined 
to  American  Gear  Manufacturers  Association  (AGMA)  class  13 
accuracy.  The  gear  profiles  were  modified  with  linear  tip  relief  as 
shown  in  Fig.  2.  Test  gear  parameters  are  shown  in  Table  I. 

Tooth  root  fillet  strains  were  measured  on  the  tensile  and  com¬ 
pressive  sides  of  two  successive  teeth,  at  the  30*  tangency  location 
(Cornell,  1980),  Fig.  3.  Dynamic  strains  were  recorded  for  the  four 
gages  at  a  matrix  of  28  load-speed  test  conditions:  four  speeds 
(800,  2000,  4000,  and  6000  rpm)  and  seven  torque  levels  (16,  31,  47. 


63,  79,  94,  and  110  percent  of  the  nominal  torque  of  71.8  N-m 
(635  lb-in.)).  The  strain  gage  signals  were  transmitted  through  a 
slip-ring  assembly  to  constant-current  signal  conditioners.  The 
dynamic  strain  data  were  digitized  by  a  digital  data  acquisition 
system,  with  sample  rates  of  50  kHz  per  channel  at  the  2000-,  4000-, 
and  6000-rpm  speeds  and  20  kHz  at  the  800-rpm  speed.  A  once-per- 
revolution  timing  pulse  provided  an  accurate  singular  position  refer¬ 
ence  for  digital  resampling  by  linear  interpolation.  The  resampled 
data  were  then  synchronously  averaged  over  several  successive  revo¬ 
lutions  to  reduce  noise  effects,  such  as  the  torque  fluctuation  from 
the  drive  belt. 

The  averaged  strain  data  values  were  converted  to  normal  tooth 
force  (dynamic  tooth  load)  by  using  calibration  data  measured 
under  static  conditions.  The  calibration  apparatus  and  data  reduc¬ 
tion  procedures  are  more  fully  described  in  Rebbechi  et  al.  (1991) 
and  Oswald  et  al.  (1991). 

ANALYTICAL  MODEL 

ARL_DYN  (Rebbechi,  1991),  Fig.  4,  considers  both  torsional 
and  lateral  displacements  of  the  gears  and  can  accommodate  vari¬ 
able  numbers  of  teeth  in  contact.  The  equations  of  motion  were 
derived  by  considering  dynamic  equilibrium  for  moments  and  forces 


(b)  Detail  of  gearbox. 
Figure  1 . — NASA  gear  noise  rig. 


acting  on  the  components  of  the  system.  Moments  are  taken  about 
the  axes  of  rotation.  The  mode  of  deflection  of  the  gear  teeth  is 
taken  as  rotation  about  a  point  in  the  gear  wheel  at  one  tooth 
height  below  the  base  circle  (see  Rebbechi,  1983).  The  tooth  stiff¬ 
ness,  which  varies  with  load  position,  is  calculated  according  to  a 
deflection  equation  in  Merritt  (1971).  Forces  acting  through  gear 
centers  are  resolved  into  coordinates  xt,  y (,  x2,  y2.  This  results  in 
[8+2(t-l)j  equations  of  motion,  where  t  is  the  number  of  tooth 
pairs  in  contact.  The  set  of  equations  can  be  reduced  to  (7-t-(t-l}| 
equations  by  invoking  the  kinematic  constraint  equations  for  the 
gears  in  mesh. 

To  improve  numerical  accuracy,  the  equations  involving  gear 
body  rotation  91 ,02  were  transformed  to  new  coordinates  ©  (  =  (flj 
+  nS2)/ 2  (average  rotation  of  gear  pair)  and  ©2  =  (0;  -  n0}  —  c)/2 
(relative  rotation  of  gear  pair),  where  61  and  02  represent  the  rota¬ 
tion  of  the  gear  wheels,  n  is  the  gear  ratio,  and  c  is  a  constant  used 
to  make  ©2  =  0  if  there  are  no  errors.  The  new  equations  relate  to 
(1)  the  absolute  rotation,  and  (2)  the  relative  rotation  of  the  gear 
wheels.  The  gear  dynamics  code  assembles  the  equations  of  motion 
in  matrix  form  as  shown  in  equation  (1): 

+  [C]{X}  -  |Kj{X}  =  {F}  (1) 

Modification, 


Roll  angle,  deg 
(b)  Driven  gear. 

Figure  2.— Test  gear  profile  traces. 
TABLE  I.— TEST  GEAR  PARAMETERS 


Gear  type . Standard  involute;  full-depth  tooth 

Number  of  teeth  .  28 

Module,  mm  (diametral  pitch,  in.-1)  .  3.175  (8) 

Face  width,  mm  (in.) .  6.35  (0.25) 

Pressure  angle,  deg  .  20 

Theoretical  contact  ratio  . 1.64 

Driver  modification  amount,  mm  (in.) .  0.023  (0.0009) 

Driven  modification  amount,  mm  (in.)  .  0.025  (0.0010) 

Driver  modification  start,  deg .  24 

Driven  modification  start,  deg  .  24 

Tooth-root  radius,  mm  (in.) .  1.35  (0.053) 

Gear  quality . AGMA  class  13 

Nominal  (100  percent)  torque.  N-m  (in. -lb)  ....  71.77  (635.25) 


2 


yy/yy 


fc»l  V. 


\ 


'V 


A 


\ 


//' 


/ 


rotation  of  gear  body  1 
rotation  of  gear  body  2 
rotation  of  gear  teeth 
lateral  displacement  of  gear 
wheel  1,  transverse  direction 
lateral  displacement  of  gear 
wheel  1 ,  radial  direction 
lateral  displacement  of  gear 
wheel  2,  transverse  direction 
lateral  displacement  of  gear 
wheel  2,  radial  direction 


r'  r*V\ 


V 


*'W 


C-89-1 0864 


(a)  Gage  installation. 
r-  Tooth  1 

30°  tangency  \  / 

to  root  fillet  \  '  /  r 

(location  of  /'N<TN\  ^  / 

gages  1  to  4)  --  pVAr®/ 


Tooth  2 


/ 

Gage 


tX 


Tooth  0- 


Rotation 


\  i 


I  \ 

(b)  Gage  location. 

Figure  3. — Strain  gage  installation  and  location  on  test  gear. 

*h«re  M,  C,  and  K  are  square  coefficient  matrices  representing  the 
mass  (inertia),  damping,  and  stiffnesses,  and  X  and  F  are  vectors  of 
displacements  and  forces  (torques),  respectively.  The  order  (number 
°f  degrees  of  freedom)  of  equation  (1)  is  (6+t),  where  t  if  the  num- 
***'  of  tooth  pairs  in  contact.  For  low-contact-ratio  gears,  C  nor¬ 
mally  varies  between  1  and  2.  If  tooth  separation  (tooth  bounce) 
^turs,  C  may  also  be  zero.  Therefore,  the  degrees  of  freedom  vary 


Figure  4. — Gear  dynamics  model. 

between  6,  7,  and  8,  depending  on  whether  there  are  zero,  one,  or 
two  pairs  of  teeth  in  contact.  The  M,  C,  and  K  matrices  contain 
a  number  of  nonlinear  terms  owing  to  the  nonlinear  kinematic 
constraints. 

The  dynamics  code  allows  prescription  of  such  features  as  profile 
errors  and  modification,  shaft  deflection  (including  interaction 
effects  with  conditions  of  tooth  contact),  tooth  deflection  (including 
resulting  change  of  contact  position  and  common  normal  direction), 
and  tooth  sliding  friction.  The  friction  coefficient  for  the  gear  mesh 
(the  friction  force  divided  by  the  normal  force)  was  taken  to  be  0.06. 
Values  up  to  0.10  were  tried  but  made  little  difference  in  the  results. 
Material  damping  in  the  gear  tooth  was  modeled  as  viscous  damping 
and  expressed  as  the  damping  ratio  (fraction  of  critical  damping). 
Gear  tooth  material  damping  produces  a  significant  effect.  A 
damping  ratio  of  0.10  gave  the  best  correlation  with  experimental 
data.  The  same  damping  ratio  (0.10)  was  assumed  for  lateral 
bending  of  the  shaft.  The  torsional  shaft  damping  coefficient  was 
2.3xl0-i  N-m/(rad/s ec).  Changing  the  shaft  damping  values  had 
little  effect. 

For  input  and  output  boundary  conditions,  steady  externa! 
torques  were  assumed.  The  code  solves  the  equations  of  motion 
with  a  Newmark-Beta  nur-ierical  integration  lecm.iquc.  The  profile 
modification  measured  for  the  test  gears  (Fig.  2)  was  specified  for 
the  analysis. 


3 


results  and  discussion 


Dynamic  tooth  loads  computed  from  the  strain  gage  readings 
were  compared  with  the  predictions  of  ARL_DYN  for  dynamic 
tooth  force.  Initial  runs  of  the  computer  model,  although  producing 
good  results  at  high  loads  (with  an  average  error  in  maximum  tooth 
load  of  5  percent),  did  not  successfully  characterise  dynamic 
behaviour  at  light  loads  and  high  speeds.  An  example  is  shown  in 
Fig.  5,  where  at  the  roll  angle  of  21°,  the  predicted  value  is  about 
3.5  times  the  measured  value. 


Figure  5. — Comparison  between  initial  prediction 
of  tooth  force  and  measured  result  at  4000  rpm 
and  31 -percent  torque. 


single-tooth  contact  to  the  tooth  tip.  It  appears  that  the  overshoot  f 
in  response  evident  in  Fig.  6  arose  from  the  inability  of  the  gear  pair  I 
to  instantaneously  adapt  to  the  “separating  mode”  (where  the  | 

driven  gear  leads  the  driving  gear).  As  the  gear  speed  increased,  » 

this  effect  became  more  marked.  The  dominant  factors  here  are  j 

probably  the  gear  wheel  inertia,  which  relates  to  speed  through  the  I 

separating  acceleration,  and  the  external  torque,  which  acts  to  \ 

reduce  dynamic  overshoot.  j 

The  effect  of  overshoot  or  tooth  separation  was  reproduced  in  > 
the  computer  model.  Initially,  the  value  used  for  gear  wheel  inertia  I 
did  not  include  the  inertias  of  the  gear  hub,  shaft,  spacers,  or  coup¬ 
lings.  When  these  parts  were  accounted  for,  the  inertia  value  in¬ 
creased  by  a  factor  of  4.  The  tooth  forces  for  the  “light”  and  the 
“standard”  (corrected)  gears  are  compared  in  Fig.  8.  The  measured 
tooth  force  data  are  also  shown  for  comparison.  This  result  closely 
accords  with  and  confirms  the  hypothesis  of  Munro  (1989),  who  des¬ 
cribed  the  tendency  of  gears  with  long  relief  to  separate  at  light  loads. 

The  influence  of  shaft  deflection  was  also  considered.  Owing  to 
the  construction  of  this  test  gearbox  (Fig.  1),  where  the  gears  are  j 
centrally  mounted  on  relatively  long  supporting  shafts,  it  was  at  one 
stage  thought  that  lateral  deflection  may  be  a  cause  of  the  tooth 
bounce  observed  when  light  loads  are  combined  with  high  speeds. 
Dynamic  tooth  strains  at  4000  rpm  are  compared  for  seven  torque 
levels  in  Fig.  9  The  effect  of  tooth  bounce  can  be  seen  in  the  curve 
for  16-percent  torque.  Here,  the  force  vanishes  around  the  pitch 
point,  indicating  that  the  teeth  have  lost  contact.  It  is  interesting 
to  note  that  the  tooth  bounce  shown  here  is  not  unique  to  the 
NASA  gear  noise  rig  but  is  also  seen  in  the  results  of  other- 
researchers  such  as  Tobe  et  al.  (1977).  Figure  10  compares  the  pre¬ 
dicted  dynamic  tooth  force  for  the  norma!  gears  with  a  case  where 
the  shaft  stiffness  is  increased  by  a  factor  of  4.  There  is  little 
difference  in  the  character  of  these  curves. 

The  shaft  deflection  in  the  radial  direction  (along  the  line 
joining  the  gear  centers)  for  both  normal  and  stiff  shafts  is  plotted 
in  Fig.  11.  As  expected,  the  mean  deflection  was  less  when  the  stiff¬ 
ness  was  increased,  but  surprisingly  the  dynamic  displacement 
increased.  This  increase  in  dynamic  displacement  for  a  stiffer  shaft 
was  apparently  a  resonance  effect  (note  the  phase  shift  of  approxi-  ’> 
mately  90°  between  the  curves).  This  indicates  that  shaft  flexi-  j 

bility  is  not  likely  to  be  a  contributing  factor  to  the  tooth  bounce.  i 


The  experimental  results  were  critically  examined  to  evaluate 
reasons  for  the  disagreement  at  light  loads  between  analysis  and 
experiment.  A  segment  of  these  results  is  shown  in  Fig.  6.  Here  the 
dynamic  tension  strains  at  31-percent  torque  are  plotted  for  four 
speeds  in  the  range  800  to  6000  rpm.  Each  curve  has  three  parts: 

(1)  a  region  where  the  strain  increases  rapidly  as  the  load  is  taken 
up  by  the  tooth  of  interest;  (2)  a  region  of  slowly  declining  static 
strain  (with  dynamic  effect  superimposed)  where  the  entire  load  is 
carried  by  a  single  tooth  pair;  this  region  lies  approximately 
between  points  A  and  B  in  Fig.  6;  and  (3)  a  region  where  the  load  is 
passed  to  the  following  tooth.  It  can  be  seen  that  the  load-sharing 
regions  (1  and  3),  and  thus  the  effective  contact  ratio,  were  virtually 
unaffected  by  speed.  Speed  had  little  effect  on  the  dynamic  load 
until  point  A  in  Fig.  6  was  reached.  The  higher  speeds  show  an 
overshoot.”  As  this  overshoot  changed  only  in  magnitude  as  the 
speed  increased,  and  did  not  change  in  angular  position,  it  most 
likely  resulted  from  a  predominantly  inertial  effect,  and  not  from  a 
combined  mass/stiffness  (resonance)  effect,  where  we  would  expect  a 
phase  shift. 

The  approximate  displacement  error  for  this  gear  tooth  profile 
when  lightly  loaded  is  as  described  by  Munro  (1989)  and  is  sketched 
ln  l"ig-  i.  Munro  terms  this  type  of  profile  correction  as  “long 
profile  relief,  where  the  relief  extends  from  below  the  high  point  of 


Dynamic  load  predictions  for  the  model  with  the  normal  . 

(heavier)  gear  wheels  are  compared  with  measured  values  in  Fig.  12.  ! 

Agreement  is  reasonable  except  in  Fig.  12(c),  at  6000  rpm  and  ( 

31-percent  torque.  The  prediction  shows  tooth  separation  at  this 
condition,  but  tooth  separation  was  actually  recorded  at  the  lower 
speed  of  4000  rpm  and  the  lighter  torque  of  16  percent  (see  Fig.  9).  , 

It  is  evident  from  these  results  that  further  refinement  of  the  model 
is  necessary  to  produce  consistent  results  across  the  whole  speed 
range.  It  is  probable  that  the  introduction  of  additional  torsional 
degrees  of  freedom,  representing  the  motor  and  the  dynamometer, 
would  aid  in  this  regard,  so  that  the  dynamic  load  increments  at 
different  speeds  could  be  magnified  or  reduced.  The  analysis  was 
particularly  successful  in  predicting  the  response  in  the  load-sharing 
region  (roll  angle  greater  than  23°  or  less  than  19°). 


The  peak  values  of  the  dynamic  load  from  both  measured  and 
predicted  data  are  compared  at  four  speeds  (800,  2000.  4000.  and 
6000  rpm)  for  the  highest  torque  level  (1 10  percent)  in  Fig.  13.  The 
predicted  and  measured  data  show  the  same  trend  ( i.e. .  a  minimum 
at  about  4000  rpm),  and  the  values  agre  within  an  aicruec  crp'r  of 
5  percent.  The  static  load  iine  drawn  in  Fig.  13  is  calculated  from 
the  external  applied  torque  of  110  percent.  The  resulting  force  of 
1894  N  (426  lb)  is  computed  from  the  torque  divided  by  the  base 
circle  radius. 
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Figure  6. — Tension  gage  strain  at  31 -percent  torque  and  four  speeds. 


Figure  7.— Relative  displacement  in  single-  and  double-tooth- 
contact  regions  and  effect  of  speed  on  dynamic  load. 


Figure  8. — Effect  of  gear  wheel  inertia  on  pre¬ 
dicted  dynamic  tooth  force  at  4000  rpm  and 
31  -percent  torque. 
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igure  9. — Dynamic  tooth  strains  at  4000  rpm  and  seven  torque  levels. 
Compressive  strains  are  shown  as  positive  for  comparison  with  tensile  data. 


Figure  1 1  .—Predicted  radial  shaft  deflection  at  4000  rpm  and 
31  -  percent  torque.  Comparison  of  results  for  two  shaft 
stiffnesses. 
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Figure  1 2. — Comparison  between  predicted  and  measured  dynamic  tooth  loads  for  normal  gear  wheels. 
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Figure  1 3. — Comparison  between  predicted  and  measured  peak 
dynamic  tooth  load  at  110-percent  torque  and  four  speeds. 
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CONCLUSIONS 


Experimental  data  for  gear  tooth  dynamic  load  were  compared 
with  predictions  from  the  Australian  Defence  Science  and  Tech¬ 
nology  Organisation  Aeronautical  Research  Laboratory’s  gear  dyna¬ 
mics  code.  The  effects  of  lateral  shaft  stiffness  and  gear  body  inertia 
were  examined  by  using  the  computer  model  to  improve  predictions 
of  gear  tooth  bounce  as  observed  at  light  loads  and  high  speeds. 

The  following  results  were  obtained: 

(1)  Peak  dynamic  load  predictions  agreed  with  measured  data 
within  an  average  error  of  5  percent  for  110-percent  torque  and 
speeds  ranging  between  800  and  6000  rpm. 

(2)  Tooth  separation  (or  bounce)  was  observed  in  the  experi¬ 
mental  data  for  light-load,  high-speed  operation.  The  computer 
model  predicted  tooth  separation  under  slightly  different  conditions. 
The  model  shows  that  this  phenomenon  is  primarily  dependent  on 
the  operating  conditions  of  speed  and  load  and  the  physical  param¬ 
eters  of  tooth  profile  and  gear  body  inertia.  An  increase  in  gear 
wheel  inertia  increases  the  likelihood  of  tooth  separation. 

(3)  The  analytical  model  was  successful  in  simulating  the 
degree  of  load  sharing  between  gear  teeth  in  the  multiple-tooth- 
contact  region. 
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